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We consider antiferromagnets breaking both time-reversal (O) and a primitive lattice transla- 
tional symmetry (T1/2) of a crystal but preserving the combination 5* = OTi/2- The S symmetry 
leads to a Z2 topological classification of insulators, separating the ordinary insulator phase from 
the "antiferromagnetic topological insulator" (AFTI) phase. This state is similar to the "strong" 
topological insulator with time-reversal symmetry, and shares with it such properties as a quantized 
magnetoelectric effect. However, for certain surfaces the surface states are intrinsically gapped with 



a half-quantum Hall effect [a^ 



e /{2h)], which may aid experimental confirmation of ^ = tt 



quantized magnetoelectric coupling. Step edges on such a surface support gapless, chiral quantum 
wires. In closing we discuss GdBiPt as a possible example of this topological class. 

PACS numbers: 03.65.Vf, 75.50.Ee, 73.43.-f, 73.20.At, 85.75.-d 
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I. INTRODUCTION 

Several insulating materials are now known experimen- 
tally to have metallic surface states as a result of topo- 
logical properties of the bulk electron wave functions li"— 
These "topological insulators"^ resulting from strong 
spin-orbit coupling were originally understood theoreti- 
cally by classifying single-electron states of materials with 
time-reversal invariance, building on previous work in 
the two-dimensional (2D) casei^ The novel metallic sur- 
face states of the three-dimensional (3D) topological in- 
sulators can be observed directly via angle-resolved pho- 
toemission spectroscopy (ARPES) and, in the simplest 
case, can be viewed as a reduced version of graphene with 
a single surface Dirac point, rather than two in the case 
of graphene, and a single spin state at each momentum 
rather than two. 

Our goal in this paper is to explain how three-di- 
mensional antiferromagnetic insulators with broken time- 
reversal symmetry can nevertheless have nontrivial fea- 
tures similar to that of the topological insulators. Along 
some planar surfaces they have gapless surface modes, 
while along others the surface is gapped and there is a 
nonzero magnetoelectric coupling from the intrinsic ma- 
terial; an experimental signature in the latter case is the 
existence of one-dimensional (ID) metallic states along 
step edges on the surface. We concentrate here on the 
conditions for an antiferromagnetic insulator to be in the 
topologically nontrivial class and on the measurable con- 
sequences at its surfaces. 

The time-reversal invariant topological insulators are 
described by Z2 topological invariants {i.e., there are only 
two possible values, "odd" and "even") that differ from 
the integer- valued topological invariants that underlie the 
integer quantum Hall effect (IQHE) in two-dimensional 
time-reversal-breaking systems. A simple picture of the 
state we discuss is obtained by starting from a nonmag- 
netic topological insulator on a bipartite Bravais lattice, 
then adding antiferromagnetic order that doubles the 
unit cell. One of the three-dimensional topological in- 



variants survives in this process. Note that this differs in 
several ways from the two-dimensional model introduced 
by Haldane on the honeycomb lattice, ^'^ which is classi- 
fied by the standard IQHE integer- valued topological in- 
variant (TKNN intcgcr^^ or "Chern number") and where 
the time-reversal breaking does not change the structural 
unit cell, which is on the hexagonal Bravais lattice. An- 
other case previously considered is a system that breaks 
time-reversal Q and spatial inversion H but preserves the 
combination OH (note that the Haldane model does not 
preserve this combination); here there are Z2 invariants 
m d = 1,2 for spinless systems and no topological invari- 
ants for spin 1/2 systems 

The basic idea in this paper is to classify crystals with 
broken time-reversal but with an unbroken symmetry 
of the form S — 0Ti/2, where T1/2 is a lattice translation 
symmetry of the "primitive" (structural) lattice that is 
broken by the antiferromagnetic order. Because the topo- 
logical invariant involves explicitly the lattice operation 
Ti/2, it is sensitive to how this lattice operation is modi- 
fied by a surface, as mentioned above, and even the gap- 
less surface state is not expected to be stable to disorder 
(in contrast to the conventional topological insulator). A 
macroscopic description is useful in order to understand 
the conditions for the topological antiferromagnet to be 
stable. The three-dimensional topological insulator can 
be characterized by the existence of a quantized magneto- 
electric coupling in the electromagnetic Lagrangiani^"— 
(c=l): 
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The coupling 6 is only defined modulo 27r, and ordinary 
insulators with time-reversal invariance have 9 = Q. The 
presence of either time-reversal symmetry or inversion 
symmetry is sufficient to guarantee that the other orbital 
magnetoelectric terms are absent fi^iii The product S is 
also enough to guarantee that the space-averaged 9 is 
quantized to zero or tt, since 9 is odd under S. 

The bulk value 9 = tt allows either metallic surfaces 
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or gapped surfaces, but in the gapped case there must 
be a half-integer quantum Hall effect. In the conven- 
tional topological insulators, the surfaces are intrinsically 
metallic and observation of the magnetoelectric coupling 
seems to require adding a time-reversal-breaking pertur- 
bation. In the topological antifcrromagnets, some sur- 
faces have a gapped state just from the material's own 
time-reversal-breaking, which suggests that experimental 
confirmation that 9 = tt, which has not yet occurred, may 
be easier in these materials, using the same techniques 
previously used to extract 9 in Cr203J^ Surface disor- 
der would complicate that approach but would enable 
observation of special features at step edges as discussed 
below. 

In the following section, we define the topological an- 
tiferromagnet in terms of band structure and verify the 
connection to the macroscopic description in terms of 
magnetoelectric response. Then the surface properties 
are discussed, which will likely be important for experi- 
mental detection. In closing we discuss the requirements 
for experiment and comment on the possibility that the 
antiferromagnetic Heusler alloy GdBiPt may realize this 
phaseji^ the possibility that such Heusler alloys may in- 
clude several topological insulators has recently been a 
topic of interest i^iSi 



satisfies T^^^O^) 
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II. 



Z2 TOPOLOGICAL INVARIANT 



In this section, we construct the Z2 invariant which 
distinguishes between the trivial insulator and "antifer- 
romagnetic topological insulator" (AFTI) phases. We 
consider a antiferromagnet breaking both the primitive 
lattice symmetry T'1/2 and time- reversal symmetry 0, 
but preserving the combination S — QTi/2- The unit 
cell is effectively doubled as a result and is the new 
lattice translation (which accounts for the notation). In 
the following, lattice vectors are elements of this doubled 
lattice except where otherwise specified. 

A free particle Hamiltonian takes the form T-L — 
X^keBZ ^k-^(^i' ^2,^3)^14 hi reciprocal space, where 
and ^' are fermion creation and annihilation operators; 
ki,k2,k^ € [0, 27r) are momentum coordinates defined by 
ki = k • a^; and are the lattice translation vectors. 
The eigenvectors Uk of the Hamiltonian -ff (k) are related 
to the wave functions by Bloch's theorem V'k = e^^ '^uu. 
Consequently, the Hamiltonian is not periodic in k, but 
rather satisfies i?(k-|-G) = e~'*^ ''7?(k)e**^ '", where G is 
a reciprocal lattice vector and r is the position operator 
in this context. Finally, we single out as such that 
gives a translation by —a.3. 

For spin-1/2 fermions, the time-reversal operator may 
be written as = — icr^/C in a suitable basis, where K, is 
the complex conjugation operator. In addition, Q (and 
S) also flips the sign of the momentum: k — k. The 
translation operator rx/2(k) will move the lattice by half 
a unit cell, so that its representation in reciprocal space 



ri/2(k) = ei' 
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where 1 is the identity operator on half the unit cell. 
Note that the operators Q and Tif2 commute so that 

eTi/2(k) = ri/2(-k)e. 

The combination ^k — 6Tx/2(k) is antiunitary like 
Q itself, but with an important difference: while 0^ = 
-1 for spin-1/2 particles, S'^ = S-^Sk = -e'''\ The 
Hamiltonian is invariant under the combination of time- 
reversal and translation: 



5kif(k)5k ' = Hi-k) 



(3) 



At the Brillouin zone (BZ) plane ^3 = the Hamilto- 
nian satisfies SH{ki, k2,Q)S~^ = H{—ki,—k2,Q) with 
(S'|fe3=o)^ = ^1- These properties lead to a Z2 topo- 
logical classification of this two-dimensional system, by 
analogy to the Z2 invariant in the quantum spin Hall 
(QSH) effect^ (the same invariant can be rederived in the 
Hamiltonian picture used here) ^ At the plane fcs = tt, by 
contrast, S'^ — +1 and there are no topological invariants 
associated with this plane. 

The Z2 invariant may be computed from the Berry 
connection and curvatur o^'^"^ on the k^ = plane, or in 
the presence of spatial inversion by looking at the four 
time-reversal momenta at ki,k2 G {0, tt}.^'* 

Even though the topological invariant is calculated 
from a two-dimensional slice in the Brillouin zone for a 
particular choice of unit cell, it reflects the topology of 
the three-dimensional band structure. For example, S 
symmetry gives no invariants in ID or 2D. In the Ap- 
pendix, we show that the 3D Z2 invariant is independent 
of unit cell choice. In the remainder of this section, we 
will give a more detailed picture of this topological phase. 



A. Relation to the time-reversal invariant 
topological insulator 

If we imagine the system described by a time-reversal 
breaking order parameter M {e.g., a staggered magne- 
tization), what happens when we restore time-reversal 
symmetry by letting M go to zero while maintaining the 
insulating phase (band gap)? 

To understand what happens, it is useful to recall 
briefly the classification of three-dimensional time-re- 
versal band insulators. In the Brillouin zone, there 
are six planes which satisfy time-reversal QH (k)Q~^ = 
i/(— k), and each has a corresponding Z2 invariant: 
ao, a^, /3o, /Stt, 7o, Ttt classify the planes ki — 0,7r, k2 — 
0, TT, and fca = 0, TT respectively. Here, we use the conven- 
tion (even) and 1 (odd) to denote the elements of Z2. 
The six values must satisfy the constraint s = ao + Q^ir — 
1^0 + Ptt = 7o + Ttt, all modulo 2, so only four combi- 
nations of these quantities are independent: s,ao,/3oi7o- 
The value s is the "strong" topological invariant, and the 
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other three Z2 are known as the "weak" invariants; to- 
gether they classify the 3D system. A strong topological 
insulator (STI) is one in which s is non-trivial, that is, 
s = 1. 

Upon doubling the unit cell in the 33 direction, the 
Brillouin zone halves by folding in the fcs direction. (In 
this subsection only, is the lattice vector of the struc- 
tural lattice, and af = 2a^ is the lattice vector of the 
"doubled" system which supports an antiferromagnetic 
coupling.) We can write the Hamiltonian of the doubled 
system in terms of the undoubled Hamiltonian i/k: 



U 



H{ki/2) 









H{ki/2 + TT) 



V2 



(4) 



Here Gg is the reciprocal lattice vector dual to af, r is 
the position operator, and the dependence on ki and k2 
are omitted for brevity. The unitary transformation U 
ensures that the eigenvectors of i?'' satisfy Bloch's theo- 
rem. 

Under the doubling process, the k^ — Q and fcg = tt 
planes collapse onto the fcg = plane. The new invari- 
ant 7q is given as a sum 70 + T-n- = s since the unitary 
transformation U does not affect any these topological 
invariants. On the other hand, the planes fca = ±7r/2 
map to the plane kf = tt. Since ±7r/2 are time-reversal 
conjugate and those planes (like all BZ planes, by as- 
sumption) have vanishing Chern numbers, it can be seen 
that 7^ is always zero. 

Adding an antiferromagnetic (0-breaking) parameter 
M to a STI produces an AFTI. Alternatively, as we turn 
down the time-reversal breaking parameter M, the anti- 
ferromagnet reverts to the doubled system. The Z2 in- 
variant describing our system is 7q = s and we have a STI 
at M = (provided the bulk gap does not close). This 
gives a natural way to construct a non-trivial topologi- 
cal antiferromagnet - by taking a STI and introducing a 
staggered magnetization which breaks time-reversal but 
preserves S. 



B. Magnetoelectric effect and the Chern-Simons 
integral 

The strong topological insulator exhibits a quantized 
magnetoelectric effect, which can be taken as its defini- 
tionJ^"— To review briefly, the magnetoelectric response 
tensor 
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is odd under the action of time-reversal. In a 0-invariant 
medium, this immediately restricts the off-diagonal ele- 
ments of the tensor to vanish. However, the ambiguity 
in defining the bulk polarization^^i2fi allows the diagonal 



elements to take a nonzero value. In fundamental units, 
the strong topological insulator has 
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with 9 ^ TT. 

The antiferromagnetic topological insulator suggested 
here does not have time-reversal symmetry microscop- 
ically; the relevant symmetry operation is S. This dis- 
tinction should not affect the macroscopic response of the 
system to uniform fields (i.e., 9), although there could be 
short- wavelength components of a* . 

From the general theory of orbital magnetoelectric re- 
sponses in band insulators, the nonzero contribution to 
a'j in cases of discrete symmetries such as time-reversal 
comes from the Chern-Simons integral. 
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where A^"^ = {uf^\id\u'Q is the Berry connection (a ma- 
trix-valued 1-form), and fi, v label filled bandsi^ The cur- 
vature 2-form is F = dA — iAl\ A. Under a gauge trans- 
formation (a unitary transformation between the bands), 
the Chern-Simons integral will change by an integer mul- 
tiple of 27r, hence only 9 mod 27r is physical. 

Under time-reversal |uk) — > ©l^k), the quantities k — > 
— k and csg — csg, and 9 changes sign. The transla- 
tion operator T^ji = e^^^^ J] changes csg by an exact 
form (total derivative) and does not affect d. Together, 
S symmetry implies that 9 = —9 ~\- ^-nn for some integer 
n, which quantizes to (topologically-trivial phase) or 
TT (topological insulator phase). 

The topological phase remains well-defined even when 
the single particle invariant is not, in the case with elec- 
tron-electron interactions. The macroscopic 9 angle re- 
mains quantized (at or tt) as long as the bulk gap does 
not close, so the AFTI is stable to sufficiently weak in- 
teractions. 

The presence of S symmetry forces the Chern numbers 
on all BZ planes to be zero. In a three-dimensional sys- 
tem, the three Chern numbers are the only obstruction 
to finding a set of continuous wave functions in the Bril- 
louin zone (respecting Bloch boundary conditions) . This 
guarantees the existence of a single- valued connection A 
for Eq. ([7]). Such A might not respect S symmetry, but 
this is no impediment to computing the Chern-Simons 
integral. 



III. CONSTRUCTION OF EFFECTIVE 
HAMILTONIAN MODELS 

In this section, we present two explicit examples of 
Hamiltonians in the antiferromagnetic topological insu- 
lator class. Henceforth, we refer these as "model A" and 
"model B." 
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A. Construction from strong topological insulators 

As noted in Sec. Ill Al we can create an antifcrroniag- 
netic topological insulator by adding a staggered time- 
reversal breaking term to a strong topological insulator. 
Here we present an explicit Hamiltonian constructed in 
such way. 

We start with a four-band model on a cubic lattice by 
Hosur et al.,^^ with four orbitals/spins per cubic site, 

H(kxi ky, kz) 

= VT^ (g) {sinkxcr^ + sin fcycr^ + sin k^a^) 

-t- [m -f t(cosfc2; -I- cos/cy -I- cosA:2)]t^, (8) 

where <t and t are two sets of Pauli matrices. This 
Hamiltonian is in the strong topological phase when 
\t\ < \m\ < 3\t\ and A 7^ 0, with the time-reversal op- 
erator represented by —ia^lC. 



To double the Hamiltonian in the z direction, first de- 
compose H{kx,ky,kz) into a hopping Hamiltonian as fol- 
lows: 

H{kx, ky, k,) = B" + Y, {B^e-'^- + B^-) , (9) 

where /i = x,y,z. The matrices describe hopping 
from adjacent cells from the direction, Bj^ are hop- 
ping from direction, and i?" describes the "self-inter- 
action" of a cell. The new lattice vectors are: 
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which defines a face-centered-cubic (FCC) lattice with 
the primitive unit cell whose volume is double that of 
the original cubic cell. 



Doubling the unit cell gives the following eight-band Hamiltonian 
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where M is a term odd under time-reversal (such as or 
T^) and represents the added antiferromagnetic coupling 
in this example. The time-reversal operator takes the 
form 



e = 









1^ (g)cr^ 



/c. 



(12) 



In the absence of M this system also has two parity 
(spatial inversion) centers, given by the operators: 



Hi = e' 









H, 









(13) 



The inversion center for Hi is between the two cubic sub- 
lattices X and Y, such that it swaps X and Y. The in- 
version center for H2 is at X, such that it takes Y to the 
next unit cell. Their product results in a translation by 
half a unit cell: H1H2 = T'i/2- 



B. Construction from magnetically induced 
spin-orbit coupling 

1. Motivation 

Consider four atoms placed in a rhombus geometry 
on the xy-plane as shown in Fig. [TJ with X and Y on 



X* 




•Y 



FIG. 1. Four atoms placed 
in a rhombus configuration 
on the a;i/-plane The cou- 
pling between X and Y de- 
pends on the magnetization 
of Ml and M2. 



opposite corners of the rhombus. In the simplest model, 
the spin-orbit coupling term from X to K is given by 
iXso X 6.2 ■ (T, where the sum is over the two paths 

X ^ Ml ^ Y, X M2 Y, and di,d2 are the 
vectors along the bonds X — )■ and — )• Y that 
the electron travels throughiiSi2S In this geometry this 
coupling vanishes as the cross products di x d2 from the 
two paths cancel. 

Now let Ml and M2 be magnetized in the +z di- 
rection. This creates a net magnetic field inside the 
rhombus breaking the symmetry between the two paths 
X M* — >■ Y . We can estimate its orbital effect by 
attaching an Aharonov-Bohm phase e^^'^ to each of the 
two paths, to produce a flux 20. The coupling from X 
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to Y now takes the form 

^^so = iAso [e"^ri x ra + e""^r2 x n] • (clfrcx) 

« 20Aso|r2 X ri|(c|,(7"cx), (14) 

where we expect to be proportional to the z-magneti- 
zation of M*. Hence the magnetization of intermediate 
sites Ml and M2 induces a spin-orbit interaction between 
the sites X and Y. 

If spins on Mi and M2 are aligned oppositely in the ±z 
direction, by contrast, there is no net magnetization in 
the rhombus and the symmetry between the two paths 
X —i' — >■ y is restored. Rotating the system by tt 
along the axis through points X and Y, taking Mi to 
M2 and vice versa, we see that there are no couplings 
between the two sites. Both cases are important in mo- 
tivating the model to follow. 



primitive cell. The electrons hop between A atoms by 
traveling through the magnetized B sites, and we can see 
that there are always two such paths A ^ B ^ A be- 
tween adjacent ^'s. From Fig. [21 it is apparent that spin- 
orbit coupling between two AVs on the same layer van- 
ishes by our argument earlier, as the intermediate sites 
have opposite magnetization. In contrast, the spin-orbit 
coupling between Al and A2 does not vanish. 



Al» Al» Al* 

1~ B t FIG. 2. Cross section of the 

model at (100) plane. The lay- 
A2' A2* A2* ^'^^ ^ unit cell are Al, B f, 

A2, B ^. Note that the magne- 
^ ^ tizations are not in-plane, but 

are only illustrated as such in 
• • Al • this figure. 



2. Effective Hamiltonian 

We start with a rock-salt (FCC) structure with the 
conventional cubic unit cell of side length 1. In this 
setup there are four 'A' sites located at (0, 0, 0) and per- 
mutations of (i, |,0), while the 'B' sites are located at 
(i, i, i) and permutations of (0,0, i). The B sites de- 
velop antiferromagnetic order along the (111) planes and 
magnetization in ±(1,1,1) direction. In the antiferro- 
magnetic state, the unit cell consists of four layers: Al, 
St, A2, and Bi. 

In this model there are spin up and spin down degrees 
of freedom at Al and A2, but the electronic degrees of 
freedom at B are eliminated, giving four "orbitals" per 



Now we describe our model with the following hop- 
ping terms: (1) spin-independent hoppings between Al 
and A2 atoms with coefhcient t, (2) spin-independent 
hoppings between A atoms on the same plane with coef- 
ficient t', and (3) spin-orbit term between Al and A2 with 
effective coupling ±A. As mentioned earlier, we take the 
energy to reside on B sites as far above the energy scales 
X,t,t', effectively eliminating those degrees of freedom in 
our model. 

We choose the primitive lattice vectors 



(15) 
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in terms of the simple cubic basis a^;, aj,, a^. The atoms Al and A2 are placed at — ^ and 
the unit cell. Written in the basis \A\), \A\), jAj), the Hamiltonian takes the form 



respectively within 



H = 



T' f/f 
U T 



T' = 2t'[cos(fci) -|-cos(fc2) +cos(fci - fc2)]l, 
where U = 2<[cos(^) + cos(fci + ^) + cos(fc2 + f-)] 

- 2iX[sin{^)<7'- + sin(fci + ^)<j^ + sin(fc2 + 



(16) 



which is gapped (in the bulk) when \t'\ < \t\,^\X\. The 
time-reversal operator has the representation 



crW 

ay 



K. 



(17) 



We are interested in the regime where t' is much smaller 
than t and A, as this leads to a gap in the surface spec- 
trum also. Unfortunately, we cannot provide a good ar- 
gument why t' (in-plane hopping) should be much less 
than t (interplane hopping) in a real material. 

This model has spatial inversion symmetry, given by 
the operator H = [J S] , which in effect swaps the lay- 



ers Al and A2. The filled bands at the momenta 
{ki,k2, h) = (0, 0, 0), (0, TT, 0), (tt, 0, 0), (tt, tt, 0) have par- 
ity — 1,— 1,— 1,+1 respectively, so the model is in the 
non-trivial topological phase. 



In this model, A is related to the parameter breaking 
time-reversal symmetry, at the same time protecting the 
bulk gap. If we turn the parameter A down to zero, we 
will not get a STI at A = 0, rather, the model becomes 
conducting. 



IV. SURFACE BAND STRUCTURE 

The bulk electronic band structure of an AFTI must be 
gapped to allow the topological distinction between the 
trivial phase and the non-trivial phase. At the boundary 
between domains of two topologically distinct phases we 
typically expect a gapless surface spectrum, as is the case 
at the edges of quantum Hall and quantum spin Hall sys- 
tems, as well as at the surfaces of the STI (vacuum is in 
the trivial phase). However, it should be noted that this 
is not strictly necessary. For example, while time-reversal 
symmetry requires doubly degenerate states, leading to 
gapless boundary modes between topological phases, it 
is known that breaking time-reversal but preserving in- 
version can give a topological phase whose surface states 
are gapped 3 

We distinguish between two classes of surfaces, de- 
pending on the plane of the cut relative to the crys- 
tal structure. We call a surface type F(erromagnetic) 
if it breaks the S symmetry in the bulk, and 
type A(ntiferromagnetic) if it preserves the symmetry. 
Heuristically, the distinction can be visualized by imag- 
ining a ferromagnetic/antiferromagnetic moment with a 
Zeeman coupling to the electron's spin, as in model A 
above. Then a type F surface will have all its spins 
aligned, and a net magnetization on the surface. A 
type A surface will have antiferromagnetic order such 
that we can always choose the primitive lattice vector 
parallel to the surface. As an example: In model A with 
staggered magnetization on a cubic lattice, {111} planes 
are type F, while planes {110} and {100} are type A. 

There are an odd number of Dirac cones on a clean 
type A surface, analogous to the STI. We can see why the 
surface (parallel to a^) is gapless by looking at the — 
line on the surface spectrum, which is the boundary of 
the ^3 = plane in the bulk BZ. Since the plane carries a 
non-trivial topological (QSH) phase, its boundary must 
be gapless. 

The Dirac cone's stability may also be explained by 
looking at a constant energy curve 7 in the surface spec- 
trum. This curve must be its own time-reversal image 
because of the symmetry between k and — k. The Berry 
phase of this curve (p = /^Tr [A] is ambiguous by integer 
multiples of 27r, so S symmetry forces this to be or tt, 
for the same reason it forced 6 = n in Sec. Ill Bl As in the 
STI, a TT phase implies that the Fermi surface encloses 
an odd number of Dirac cones. However, any defect or 
impurity will break the translational and S symmetry on 
the surface, thereby opening a gap. This is analogous to 
the effect of a magnetic defect on the surface of a STI. 

For a type F surface, S symmetry is broken on the 
surface and the usual protection for Dirac cones or con- 
ducting surfaces no longer exists. If the bulk and surface 
spectrum are fully gapped (i.e., not a semi- metal), then 
the surface will exhibit the half-integer quantum Hall ef- 
fect, to be discussed in the next section. 

In Fig. 131 we present the band structure of model A for 
slabs with type A and type F surfaces. Since this model 



6 




FIG. 3. (Color online) The bulk and surface band structure 
for model A, along the (100) (type A, parallel to 3.2, as) and 
(III) plane (type F, parallel to ai,a2) respectively. The red 
dots indicate surface modes. The parameters used are: v = 
0.5, m = 2, f = 1, M = cr^ with 13 layers. 



is built from a STI, the band structures are similari^ 

For model B, the surface parallel to ai and a2 is type F 
and its excitations are exactly solvable, with dispersion 

Es(ki,k2) = 2t'(cosfci + cosfc2 -I- cos(fci - ^2)). (18) 

As the dispersion shows, a surface spectrum exists for 
all values of t' and non-zero values of t and A. This 
model has the peculiar feature that the surface spectrum 
is completely disconnected from the bulk, that is, it forms 
a complete two-dimensional band structure. Figure |4] 
shows the bulk and surface band structure for two dif- 
ferent cuts. In the (111) cut, a small t' is desired if we 
want to avoid band overlaps between the valence, con- 
ductance, and surface spectrum, giving us an insulator. 

In the presence of a sufficient number of random de- 
fects, we expect that the surface electronic states are de- 
scribed by the unitary symmetry class because of the 
broken time-reversal symmetry. That symmetry class 
only has extended states at isolated values of the chem- 
ical potential; in general the surface state will have zero 
diagonal conductivity, with half-integer quantum Hall 
plateaus. The transitions between these plateaus appear 
when the chemical potential passes through an extended 
state. These transitions can be regarded as a realiza- 
tion of the two-dimensional quantum Hall effect in zero 
net field discussed by Haldane.^- Note that since both 
top and bottom surfaces of a slab will have half-integer 
plateaus, the total quantum Hall effect when diagonal 
conductivity is zero is always integral, as required for a 
single-electron two-dimensional system. 
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FIG. 4. (Color online) The bulk and surface band structure 
for model B, along the (111) plane (type A, parallel to a2, ag) 
and (111) plane (type F, parallel to ai,a2) respectively. The 
red dots indicate surface modes. The parameters used are: A = 
0.5, t=l,t' = 0.1 with 13 layers. 



V. FERROMAGNETIC SURFACES AND 
HALF-INTEGER QUANTUM HALL EFFECT 

In this section we present two perspectives on the half 
quantum Hall effect on type F surfaces, along with nu- 
merical calculations to justify our claim. 

If one views the antiferromagnet as a STI with time- 
reversal breaking term opening a surface gap, then the 
half QHE can be viewed as the root of the bulk magne- 
toelectric coupling 9 = tt. This effect follows from the 
gapped Dirac dispersion of the surface states. The sign 
of the Hall conductance depends on the sign of the effec- 
tive Dirac masS )^°'^^i'^^ which here is set by the direction 
of the Zeeman field at the surface. 

An alternate perspective of the AFTI surface comes 
from a comparison to the quantum spin Hall effect. As 
described in Sec. |lTl the Z2 invariant is computed from 
the two-dimensional plane fca = and the symmetry op- 
erator 5'|fc3=o in precisely the way that the quantum spin 
Hall (QSH) invariant is computed from the two-dimen- 
sional BZ and 8. The QSH insulator may be constructed 
by combining two copies of a QH layer with opposite 
spin and Chern number ±n. Time-reversal takes one 
layer to the other, making the combination of the two Q- 
invariant. In reality spin is rarely conserved, allowing the 
two layers to mix, making the Chern number of each spin 
ill-defined. However, a residual Z2 topological invariant 
remainS)2i2^ and we can consider the QSH as being adia- 
batically connected to the two-QH-layer model, but with 
the topological invariant n only preserved mod2. 

By analogy, we can construct an AFTI by stacking QH 
layers, with alternating Chern number of ±1 (Hall con- 
ductivity ie^/Zi), as shown in Fig. [5] The "up" (-fl) 



FIG. 5. Construction of antifer- 

I 1 romagnetic topological insulator 

3-3 r ' —I by staggering quantum Hall lay- 

/ !***** _! ers. The shaded and unshaded 
I « » » » » I boxes represents Chern number 

I . « . ,~r| of ±1. The left and right sur- 

I 1 faces are type A and gapless, 

I I while the top and bottom sur- 

— - — - — - — faces are type F with half-quan- 
tum Hall effect. 



layers are related to the "down" (—1) layers by S sym- 
metry, hence they are spatially offset from one another. 
Just like the QSH case, we can expect the layers to couple 
to one another, in a way that makes the Chern number 
ill-defined on a per-layer basis. Once again, it is appro- 
priate to consider the AFTI to be adiabatically connected 
the staggered QH layer model. In the stacked QH model, 
the Hall conductance in the bulk averages to zero, as the 
conductance of any individual layer is cancelled by neigh- 
boring layers of opposite type. In other words, any long- 
wavelength probe of the system will be unable to discern 
the individual QH layers. However, the QH layers at ei- 
ther end of the stack are not completely cancelled, there 
is a half QHE at both surface. 

To confirm this picture, we can consider a slab with 
type F surfaces and compute the 2D Hall conductivity as 
a function of position (layer). In units of e^//i, the (two- 
dimensional) conductivity in layer n can be computed 
from^^ 



2-K 



tr 



V{dV)hVn{dV) 



(19) 



Here V = X^occ l"k)(Mk| is the projector onto occupied 
wave functions at k and Vn is the projector onto basis 
states localized in layer n. 

C(n ) 





5 10 15 20 25 3 



FIG. 6. (Color online) Hall conductance spatially resolved for 
32 (blue) and 33 (red) layers of model B. The Fermi level is 
set below zero to include only the bulk valence modes and 
no surface modes. The parameters used for this plot are: 
A = 0.5, t = l,t' = 0.1. 

Figure [S] shows the results of such a computation on 
a slab cut from the rock-salt model (B) introduced in 
Eq. with type F surfaces. In this model, when the 
Zeeman field on opposite surfaces points in opposite di- 
rections (blue, upper curve) the total conductance of the 
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slab is C = 1, with each surface having a net C = 1/2; 
adding a layer such that the two surfaces have the same 
Zeeman field switches the conductance on that surface 
from +1/2 to —1/2, so that the total slab conductance 
vanishes. Note that the total conductance of a slab is 
always an integer, as requiredj^^^i^ 

Now, at the interface between two integer quantum 
Hall domains whose conductance C differs by 1, there 
will be a chiral boundary mode with conductance e^//i, 
which can be thought of as "half a quantum wire." In 
the situation outlined above, putting the two slabs with 
different conductance together is equivalent to making a 
slab with a step edge on one surface, and the chiral mode 
will reside at this step edge. Such a mode should give an 
observable signature in a tunneling experiment (FiglT]). 




FIG. 7. (Color online) ID quantum wire on type F surface 
step edge. The red and blue regions represent ferromagnetic 
layers magnetized in opposite directions. There is a gapless 
chiral quantum wire at each step edge, with chirality indicated 
by the arrow at the edge. 

It is natural to ask, what if one rotates the antiferro- 
magnetic moment by tt, flipping all the spins and effec- 
tively "peeling" off a layer of type F surface? Since the 
sign of the surface conductance C changes during this 
process, the surface (or bulk) gap must close at some 
magnetization orientation. This is analogous to applying 
a magnetic field to a STI surface. For B parallel to the 
surface, the Dirac cone shifts in momentum space but 
no gap opens. Any infinitesimal component of B out of 
the plane will open a gap, hence going from B out of the 
surface to B into the surface must necessarily close the 
surface gap. (In model B, the bulk gap would close while 
rotating the magnetization.) 

VI. CONCLUSIONS AND POSSIBLE 
RELEVANCE TO GdBiPt 

In this paper, we have looked at the topological clas- 
sification of materials breaking both time-reversal 9 and 
translational symmetry T1/2, but preserving the combi- 
nation S = QTi/2, and found a Z2 classification within 
the S symmetry class that leads to the existence of an an- 
tiferromagnetic topological insulator (AFTI) . In the most 



basic picture, an AFTI can be obtained from adding a 
staggered magnetization to a strong topological insulator 
(STI). Macroscopically, S symmetry implies a quantized 

magnetoelectric response f5|5_Q = with 9 = tt 

for an AFTI. We have also demonstrated that the sur- 
face spectrum depends on the surface cut, classified as 
type A/F. Type A surfaces possess an antiferromagnetic 
order that preserves S symmetry, with associated gapless 
excitations that can be gapped by disorder. Type F sur- 
faces break S symmetry and are typically gapped, anal- 
ogously to the situation of a Zeeman field on the surface 
of a STI. The new AFTI state is topological in a weaker 
sense than the strong 3D topological insulator, because 
its surface state is dependent on the surface plane and 
not generally stable to disorder; in that respect it is sim- 
ilar to the weak topological insulator in 3D or the "Hopf 
insulator."— (The number of Dirac cones in a STI also 
depends on the surface plane, but there is always an odd 
number of such cones.) 

The magnetoelectric coupling ^ n requires the half 
quantum Hall effect at the surface, provided the surface 
spectrum is gapped. Our numerical calculations based 
on explicit band models agree with these results. Finally, 
we predict the existence of chiral ID quantum wires at 
type F surface step boundaries, an experimental signa- 
ture verifiable via scanning tunneling measurements. 

The recent proposals that many Heusler compounds 
may be topological insulators j^i^ together with the an- 
tiferromagnetic order in GdBiPt below 9 suggest a 
possible candidate for the state proposed here. Trans- 
port experiments indicate that GdBiPt is a semiconduc- 
tor with a narrow gapJ^^ The Gd sites form an FCC lat- 
tice and hence their antiferromagnetic interaction is frus- 
trated, and further experiments {e.g. neutron scattering) 
are required to determine if the antiferromagnetic order 
falls under the S symmetry class described in this paper. 
At least one related Heusler antiferromagnet (MnSbCu) 
is known to have antiferromagnetic ordering [alternating 
spin directions on (111) planes] which belongs in the S 
symmetry class.— Should the material be truly insulat- 
ing (i.e., have a bulk gap) through its antiferromagnetic 
transition, it suffices in principle to check if it is a strong 
topological insulator above the Neel temperature. 

We have provided a topological classification and ex- 
perimental consequences for a particular combination of 
time-reversal symmetry and a lattice symmetry (0Ti/2). 
Other such combinations of time-reversal and crystal 
symmetries could lead to new topological materials be- 
yond those in the exhaustive classification of topological 
insulators stable to disorder 



VII. ACKNOWLEDGEMENTS 

The authors gratefully acknowledge discussions with 
S. Ryu, O. Yazyev, and D. Vanderbilt. The work was 
supported by NSF under Grant No. DMR-0804413 (R.M. 
and J.E.M.) and WIN (A.E.). 



9 




y — y>A 




9 









9 





FIG. 8. (Color online) Example of changing the unit cell, 
(a) The original cell (rectangle with dashed border) is trans- 
formed to the new cell (shaded rectangle with solid bor- 
der) by keeping the X portion fixed and changing Y , where 
R = — a2 — as is the displacement vector. In this example, the 
original vectors a2 = —x ~ y, as — 2y, and ai points out of 
the plane. The new vectors as = as -|- 2R = — 2a2 — a.^ = 2x 
and we choose ai and a2 to remain fixed, (b) The Brillouin 
zone. The blue plane is the effective Brillouin zone (EBZ) 
for fcs = 0, and the red plane is the EBZ for fca = 0. The 
Z2 invariant computed for these two planes are the same, (c) 
Deformation of the new EBZ (red), which decomposes into 
the old EBZ (blue) and a boundaryless plane (green). 



Appendix A: Invariance to choice of unit cell 

The construction of the Z2 invariant in Sec. HIl required 
a certain choice of unit cell. In this section, we will 
demonstrate that different choices of the unit cell will 
yield the same result. In particular, we show that differ- 
ent ways to choose the doubled unit cell are equivalent 
given a choice of structural cell. 

Begin with a Hamiltonian H defined for a set of primi- 
tive translation vectors a^, along with the operators T1/2 
such that translate by —as. We can always divide 
the Hilbert space in to two subspaces: X and Y, such 
that the translation operator Tif2 takes Y to X, and X 
to the Y in another unit cell. Physically X and Y rep- 
resent the structural unit cell whose symmetry is broken 
by antifcrromagnctism. 

Construct a new unit cell by leaving X fixed but taking 
Y from a cell R relative to the original. In the new 
system the lattice vector 3.3 = + 2R such that 2^/2 
translates the system by —a^. We want to show that the 
Z2 invariant calculated for the new Hamiltonian {H on 
the fcs = plane) is identical to that of the original one 
{H on the = plane). 



Here we remind the reader of the method used in this 
section to compute the Z2 topological invariant. ^'^"^ First 
we pick an "effective Brillouin zone" (EBZ) which is half 
of the Brillouin zone such that time-reversal will map it 
to the other half. The boundary of EBZ must be time- 
reversal image of itself. The element of Z2 is computed 
by the integrating the connection and curvature: 



D 



1 

2^ 



A- 

9EBZ 



EBZ 



mod 2, 



(AI) 



where the [t/(l)] connection A = X^occ {'u-\'id\u) is 
summed over occupied bands and curvature = dA is its 
exterior derivative in momentum space. The curvature 
T is "gauge invariant" (does not depend on the choice 
of basis functions for occupied states), but A depends on 
a particular choice of gauge for the wave functions. The 
boundary integral in the formula above requires that the 
wave functions at k and — k be S'-conjugate pairs. Any 
choice of the EBZ will give the same Z2 invariant. 

The effect of the coordinate transformation 
^1,^2,^3 — ^ ^1,^2,^3 changes the EBZ on which 
we compute the topological invariant. Since the momen- 
tum variables are related by — -\- 2R • k, we can 
always choose the EBZ for the new and old systems such 
that they share a common boundary, namely, the two 
lines satisfying R • k e {0, tt}. This guarantees that the 
boundary integral terms {§A) in (Eq. lAip are identical 
in the two cases. 

As for the term integrating curvature over the EBZ, we 
can try to deform the new EBZ to match the old EBZ. 
This deformation is allowed by the fact that J-' — dA 
is a closed 2-form; any local deformation to the surface 
{i.e., one that preserves A on the boundary) will preserve 
the integral / J-. As Fig. 8(c) shows, we cannot always 



deform one EBZ to the other; however, we can always 
decompose the new EBZ into the old EBZ plus planes 
with no boundaries. These closed planes which are either 
contractible, or they span a torus in the Brillouin zone. 
S symmetry requires that the Chern number vanishes on 
all closed two-dimensional surfaces, and it follows that 
the integral (Eq. lAll) evaluates to the same value for new 
and old unit cell. In other words, the Z2 invariant does 
not depend on how we choose the unit cell. 

We can also view the Z2 invariant as an obstruction 
to finding a continuous basis (along with the appropri- 
ate Bloch periodic boundary conditions) for the wave 
functions respecting S symmetry in the entire Brillouin 
zone.^'^ The material is in a trivial phase if such a basis 
exists. This intepretation is much harder to "compute" 
then the original definition, but is powerful in what it 
implies. For example, any (single-valued) unitary trans- 
formation or a change of coordinates will not affect the 
obstruction of finding such basis, and it is rather straight- 
forward from the definition that the Z2 invariant is inde- 
pendent of unit cell choice. 
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